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Abstract 

We investigate matter-induced (or extrinsic) CPT violation effects in neutrino oscillations in 
matter. Especially, we present approximate analytical formulas for the CPT-violating probability 
differences for three flavor neutrino oscillations in matter with an arbitrary matter density profile. 
Note that we assume that the CPT invariance theorem holds, which means that the CPT violation 
effects arise entirely because of the presence of matter. As special cases of matter density profiles, 
we consider constant and step-function matter density profiles, which are relevant for neutrino 
oscillation physics in accelerator and reactor long baseline experiments as well as neutrino factories. 
Finally, the implications of extrinsic CPT violation on neutrino oscillations in matter for several 
past, present, and future long baseline experiments are estimated. 
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I. INTRODUCTION 
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have been performed in order to incorporate the so-called LSND anomaly 



|22| within the description of standard three flavor neutrino oscillations. However, 



ihis req uires a new mass s quar ed difference different from the ones coming from atmospheric 
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one would need to have three mass squared differences instead of two - a scenario, which is 
not consistent with ordinary models of three flavor neutrino oscillations. Therefore, in most 



of the studies on CPT violation 



15 



different mass squared 
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differences and mixing parameters are introduced phenomenologically by hand for neutrinos 
and antineutrinos. This results, in the three neutrino flavor picture, in two mass squared 
differences and four mixing parameters for neutrinos and the same for antineutrinos, i.e., 
in total, four mass squared differences and eight mixing parameters. Thus, it is possible 
to have a different mass squared difference describing the results of the LSND experiment 
other than the ones describing atmospheric and solar neutrino data. It should be noted that 
;he results of the LSND experiment will be further tested by the MiniBooNE experiment 
3(|, which started running in September 2002. Furthermore, it should be mentioned that 
the standard way of incorporating the LSND data is to introduce sterile neutrinos, and 
therefore, the introduction of fundamental CPT violation, sometimes also called genuine 
CPT violation, serves as an alternative description to sterile neutrinos. However, neutrino 
oscillations between pure sterile flavors and active and sterile flavors have, in principle, been 
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401 ]. a milestone of local 



excluded by the SNO experiment 

In CPT violation studies, the CPT invariance theorem 3^, 
quantum field theory, obviously does not hold, and in addition, fundamental properties such 
as Lorentz invariance and locality may also be violated. However, the SU(3) x SU(2) x U(l) 
Standard Model (SM) of elementary particle physics, for which the CPT theorem is valid, is 
in very good agreement with all existing experimental data. Therefore, fundamental CPT 
violation is connected to physics beyond the SM such as string theory or models including 
extra dimensions, in which CPT invariance could be violated. 

The recent and the first results of the KamLAND experiment 41], which is a reactor 



long baseline neutrino oscillation experiment measuring the u e flux from distant nuclear 
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reactors in Japan and South Korea, strongly favor the large mixing angle (LMA) solution 
region for solar neutrino oscillations and the solar neutrino problem |42|. Therefore, they 
indicate that there is no need for fundamental CPT violation, i.e., having different mass 
squared differences for solar neutrinos and reactor antineutrinos. Thus, solar neutrino data 
and KamLAND data can be simultaneously and consistently accommodated with the same 
mass squared difference. 

In this paper, we investigate matter-induced (or extrinsic) CPT violation effects in neu- 
trino oscillations in matter. In a previous paper |43j, the interplay between fundamental and 
matter-induced T violation effects has been discussed. In the case of CPT violation effects, 
there exists no fundamental (or intrinsic) CPT violation effects if we assume that the CPT 
theorem holds. This means that the matter-induced CPT violation is a pure effect of the 
simple fact that ordinary matter consists of unequal numbers of particles and antiparticles. 
Matter-induced CPT violation, sometimes also called fake CPT violation, have been stud- 
led and iUustrated in some papers Q H Q H Q 0, in which numerical calculations of 
CPT-violating asymmetries between survival probabilities for neutrinos and antineutrinos 
in different scenarios of atmospheric and long-baseline neutrino oscillation experiments have 
been presented. Here we will try to perform a much more systematic study. 

The paper is organized as follows. In Sec. m we discuss the general formalism and 
properties of CPT violation in vacuum and in matter. In particular, we derive approximate 
analytical formulas for all CPT-violating probability differences for three flavor neutrino 
oscillations in matter with an arbitrary matter density profile. The derivations are performed 
using first order perturbation theory in the small leptonic mixing angle #13 for the neutrino 
and antineutrino evolution operators as well as the fact that Am^ -C Anijj ~ Am^, i. e., the 
solar mass squared difference is some orders of magnitude smaller than the atmospheric mass 
squared difference. At the end of this section, we consider two different explicit examples 
of matter density profiles. These are constant and step-function matter density profiles. In 
both cases, we present the first order perturbation theory formulas for the CPT probability 
differences as well as the useful corresponding low-energy region formulas. Next, in Sec. 
we discuss the implications for long baseline neutrino oscillation experiments and potential 
neutrino factory setups as well as solar and atmospheric neutrinos. We illuminate the 
discussion with several tables and plots of the CPT probability differences. Then, in Sec. IIVI 
we present a summary of the obtained results as well as our conclusions. Finally, in App. El 



3 



we give details of the general analytical derivation of the evolution operators for neutrinos 
and antineutrinos. 



II. GENERAL FORMALISM AND CPT- VIOLATING PROBABILITY DIFFER- 
ENCES 

A. Neutrino Oscillation Transition Probabilities and CP, T, and CPT Violation 



Let us by P[y a — > vp) denote the transition probability from a neutrino flavor a to 
a neutrino flavor (3, and similarly, for antineutrino flavors. Then, the CP, T, and CPT 
(-violating) probability differences are given by 
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and T probabilit y d ifferences have prev iously been exten- 




. In this paper, we will study in detail the CPT probability differences. 
Let us first dicuss some general properties of the CPT probability differences. In general, 
i.e., both in vacuum and in matter, it follows from conservation of probability that 
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In words, the sum of the transition probabilities of a given neutrino (antineutrino) flavor 
into neutrinos (antineutrinos) of all possible flavors is, of course, equal to one, i.e., the 
probability is conserved. Using the definitions of the CPT probability differences, Eqs. (j3J) 
and (j3J) can be re-written as 
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Note that not all of these equations are linearly independent. For example, for three neutrino 
flavors, Eqs. © and (JJJ) can be written as the following system of equations 
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A pCPT A pCPT A pCPT = 
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A pCPT A pCPT A pCPT = 
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(13) 



Hence, there are nine CPT probability differences for neutrinos and six equations relating 
these CPT probability differences. The rank of the corresponding system matrix for the 
above system of equations is five, which means that only five of the six equations are linearly 
independent. Thus, five out of the nine CPT probability differences can be expressed in terms 
of the other four, i.e., there are, in fact, only four CPT probability differences. Choosing, 
e.g., AP^p T , AP e ^ PT , AP^ PT , and AP^ PT as the known CPT probability differences, the 
other five can be expressed as 
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Furthermore, the CPT probability differences for neutrinos are related to the ones for an- 
tineutrinos by 

AP Q CPT = P{u a -> up) - P(9f3 - Pa) 

= -{P{p p - p a ) - P(i/ a - 14,)) = -APjl PT , (19) 

where a, /3 = e, //, r, . . .. Thus, the CPT probability differences for antineutrinos do not give 
any further information. 

For completeness, we shall also briefly consider the case of two neutrino flavors. In this 

5 



case, we have 
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from which one immediately obtains 

AP e c e PT = AP 0PT = -AP° PT = -AP^ T - (24) 

Thus, for two neutrino flavors there is only one linearly independent CPT probability dif- 
ference, which we, e.g., can choose as AP 0PT . 

Generally, for the T probability differences, we have j3, 

AP P = AP P = AP P = 0, (25) 
API = API = API = -API = -API = -APj T (26) 

for three neutrino flavors and 

API = API = API = AP^ = (27) 

for two neutrino flavors. Thus, in the case of three neutrino flavors, there is only one linearly 
independent T probability difference, whereas in the case of two neutrino flavors, neutrino 
oscillations are T-invariant irrespective of whether they take place in vacuum or in matter. 

Using the definitions (fT|) - (|3*)l. one immediately observes that the CP probability differ- 
ences are directly related to the T and CPT probability differences by the following formulas 



APS? + APj = AP a CPT and AP ? + AP P = AP5f T (28) 
In vacuum, where CPT invariance holds, one has APS? T = AP9i T = 0, which means 

' ' ctp ap 

that AP^p = -AFj^ and AP5j? = — AP^. Furthermore, using again the definition (JTJ), 
one finds that APf? = -AP. ?. Thus, AP ? = AP P and AP CP = AP T ,, i.e., the 

ap a fj ) ap ap a fj a fji J 

CP probability differences for neutrinos (antineutrinos) are given by the corresponding T 
probability differences for neutrinos (antineutrinos). However, in matter, CPT invariance 
is no longer valid in general, and thus, one has AP^ T ^ 0, which means that we need to 
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know both the T and CPT probability differences in order to determine the CP probability 
differences. Moreover, in vacuum, it follows in general that P{y a —>■ up) = P{pp — > v a ) and 
in particular that P{y a — > u a ) = P{y a — > u a ), which leads to AP^f = 0. Therefore, CP 
violation effects cannot occur in disappearance channels (u a — > u a ), but only in appearance 
channels (u a — > vp, where a ^ (3) 52j, while in matter one has in general AP^ ^ 0. 

In the next subsection, we discuss the Hamiltonians and evolution operators for neutrinos 
and antineutrinos, which we will use to calculate the CPT probability differences. 

B. Hamiltonians and Evolution Operators for Neutrinos and Antineutrinos 

If neutrinos are massive and mixed, then the neutrino flavor fields v a , where a = 
e, n, r, . . ., are linear combinations of the neutrino mass eigenfields u a , where a = 1,2, 3,..., 
i.e., 

n 

v a = s ^U aa v a , a = e,fjL, t, . . . , (29) 

a=l 

where n is the number of neutrino flavors and the Una's are the matrix elements of the 



unitary leptonic mixing matrix U. 



0H 



3411 Thus, we have the following relation between the 



neutrino flavor and mass states 

y «) = ^2 U aaWa), OL = 6, fl, T, . . . , (30) 



a=l 



where v a is the ath neutrino mass state for a neutrino with definite 3-momentum p, energy 
E a = \/rn\ + p 2 ~ p + ^ (if m a p), and negative helicity. Here m a is the mass of the 
ath neutrino mass eigenstate and p = |p|. Similarly, for antineutrinos, we have 

n 

\v a ) = ^2U aa \u a ), a = e,fi,r,.... (31) 

0=1 

In the ultra-relativistic approximation, the quantum mechanical time evolution of the neu- 
trino states and the neutrino oscillations are governed by the Schrodinger equation 

i±W{t)) = je{tMt)), (32) 

where \v(t)) is the neutrino vector of state and J$?(t) is the time-dependent Hamiltonian 
of the system, which is different for neutrinos and antineutrinos and its form also depends 



7 



on in which basis it is given (see App. El for the different expressions of the Hamiltonian) . 
Hence, the neutrino evolution (i.e., the solution to the Schrodinger equation) is given by 

\v(t))=e*& je W\v{t )), (33) 

where the exponential function is time-ordered. Note that if one assumes that neutrinos are 
stable and that they are not absorbed in matter, then the Hamiltonian J$?(t) is Hermitian. 
This will be assumed throughout this paper. Furthermore, it is convenient to define the 
evolution operator (or the evolution matrix) S(t,t ) as 

\v(t)) = S(t,t )\"(to)), S(t,t ) = e"^ W, (34) 
which has the following obvious properties 

S(t,t ) = S , (t,ti)S , (ti,t ), (35) 
S(t ,t ) = 1, (36) 
S(t,t )S(t,t o y = 1. (37) 

The last property is the unitarity condition, which follows directly from the hermiticity of 
the Hamiltonian J^(t). 

Neutrinos are produced in weak interaction processes as flavor states \v a ), where a = 
e,fi,r, . . .. Between a source, the production point of neutrinos, and a detector, neutrinos 
evolve as mass eigenstates |z/ a ), where a = 1,2,3, . . ., i.e., states with definite mass. Thus, 
if at time t = to the neutrino vector of state is \v a ) = \v a (to)), then at a time t we have 

n 

K(t)) = J2i S ^ t °)laKaWa)- (38) 
a=l 

The neutrino oscillation probability amplitude from a neutrino flavor a to a neutrino 
flavor (3 is defined as 

n 

A aP = {vp\v a {t)) = ^2Up a [S(t,t )] aa U* a , a,f3 = e,fi,T,.... (39) 

a=l 

Then, the neutrino oscillation transition probability for v a — * z/g is given by 

n n 

P(u a - up) = \A aP \ 2 = Y,H U ^UpaU ab U; b [S(t,t )] aa [S(t,t Q )]* bb , (40) 

a=l b=l 

where a, (3 = e, fi,r, . . .. 



The oscillation transition probabilities for antineutrinos are obtained by making the re- 
placements U aa — > U* a and S(t,to) — > S(t,t ) [i.e., V{t) — > — V(t), where V(t) is the matter 
potential defined in App. E], which lead to 

n n 

P(p a - P/j) = Y,H U ^ U k U ^ U ^ [S(t,t Q )] aa [S(t,t )]* bb = {a^b} 

a=l 6=1 

n n 

= E E KaUpaU ab U; b [S(t, t )] 1 , to)] 66 , (41) 

tt=l 6=1 

where a, j3 = e, [i,r, . . .. 

In the next subsection, we calculate the CPT probability differences both in vacuum and 
in matter. 



C. CPT Probability Differences 

In vacuum, the matter potential is zero, i.e., V(t) =0 Vt, and therefore, the evolution 
operators for neutrinos and antineutrinos are the same, i.e., S(t,t ) = S(t,t ) = e -ljHmL , 
where H m = diag {E%, E2, ■ ■ ■ , E n ) is the free Hamiltonian and L ~ t — to is the baseline 
length. Note that the Hamiltonians in vacuum for neutrinos and antineutrinos are the same, 
since we have assumed the CPT theorem. Thus, using Eqs. (J4Uj) and (J41j) . it directly follows 
that 

AP a c J T = P(u a -> up) - P(% -> u a ) = 0, (42) 

which means that there is simply no (intrinsic) CPT violation in neutrino oscillations in 
vacuum. Note that this general result holds for any number of neutrino flavors. Furthermore, 
note that even though there is no intrinsic CPT violation effects in vacuum, there could be 
intrinsic CP and T violation effects induced by a non-zero CP (or T) violation phase 5cp, 
which could, if sizeable enough, be measured by very long baseline neutrino oscillation 
experiments in the future 96]. 



In matter, the situation is slightly more complicated than in vacuum. However, the 
technique is the same, i.e., the extrinsic CPT probability differences are given by differences 
of different matrix elements of the evolution operators for neutrinos and antineutrinos. 

The probability amplitude of neutrino flavor transitions are the matrix elements of the 
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evolution operators: 

A [y a - up) = [S(t, t % a = [S f (t, t )] pa , (43) 
A(u a - v p ) = [S(t,t o )] 0a = [Sf(t,t )]p a - (44) 

Thus, we have the extrinsic CPT probability differences 



\pCPT 



[Sf(t,t )]p a [Bf(t,toj\af> ■ ( 45 ) 



In the case of three neutrino flavors with the evolution operators for neutrinos and an- 
tineutrinos as in Eqs. (jA35|) and ([ A40j) . respectively, the different AP^ T 's are now easily 
found, but the expressions are quite unwieldy. The CPT probability difference AP e ^ PT to 
first order in perturbation theory is found to be given by (see App. [X] for definitions of 
different quantities) 

a nCPT ^ |c |2 10 |2 I 1 2 I — 1 2 |5|2 

2 ^ Sil1 
= COS \ I + 



(cos29 12 5(t - t ) - jf V(f )d*'J 

- cos 2 Q - [cos29 l2 5{t - t ) + J V{t')dt'J 

1 /sin 2 f2 sin 2 f2\ . „ „ r9 . , 9 . . 

= — Jsin 2 2^ 2 (i-t ) 2 , (46) 

which is equal to zero in vacuum, in which V(t) — Vi. Note that in the case of T violation 
all diagonal elements, i.e., AP P Q , where a = e, //, r, are trivially equal to zero [c/., Eq. (J2*51) ]. 
This is obviously not the case for CPT violation if matter is present. Similarly, we find 

AP e c / T ^ \S ffil \ 2 -\S f ^\ 2 =\c w p*+is 23 fC\ 2 -\c w p-is w fA\ 2 

= c 2 3 (i/5i 2 -i/5i 2 )+4 (i^r-^i 2 ) 

+ is 23 c 23 (/3/C - /T/*C* + /3*fA - M*A*) , (47) 
AP e c / T ~ \S m \ 2 - \S L13 \ 2 = \s 23 p* - ic 23 fC\ 2 - I - s 23 (3 - ic 2 JA\ 2 
= c 2 23 {\C\ 2 -\A\ 2 ) + s 2 23 {W-\(3\ 2 ) 

- is 23 c 23 (PfC - p*f*C* + P*fA - Pf*A*) , (48) 
AP M c e PT ~ \S L12 \ 2 - \S m \ 2 = \c 23 (3 - \s 23 fA\ 2 - \c 23 (3* + \s 23 fC\ 2 
= c 2 23 (\(3\ 2 -\f3\ 2 )+4s (\A\ 2 -\C\ 2 ) 

+ is 23 c 23 (Pf*A* - 0* f A - PfC + P*f*C*) , (49) 
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AP r c e PT ~ \S f , 13 \ 2 - \S m \ 2 = | - s 23 (3 - ic 23 fA\ 2 - \s 23 p* - ic 23 fC\ 2 
= c 2 23 (\A\ 2 -\C\ 2 ) + s 2 23 (W-\(3\ 2 ) 

- ls 23 c 23 (f3f*A* - f3*fA - pfC + P*f*C*) , 

AP£f T * fe 2 | 2 - \S f , 22 \ 2 = l4«* + 4/ - i*23C 23 /(5 + D)\ 2 

- \c 2 23 a* + s 2 23 f-is 23 c 23 f(B + D)\ 2 

= 4 (M 2 - l«D - is2a4 " «T^* + " «7*£>* 

- a/P + «7*P* - afD + a*f*D*) 

+44 («/ + "7* + l^l 2 + l^l 2 + BD* + B *D 

- af-a*f* - \B\ 2 - \D\ 2 - BD* - B*D) 
-i4c 23 (B - B* + D - D* - B + B* - D + D*) , 

AP£ PT ~ \S f , 32 \ 2 - \S L23 \ 2 = | - s 23 c 23 (a* -f)+if (s 2 23 B - c 2 23 D) \ 2 

- I - s 23 c 23 (a* -f)- if (c\ 3 B - sj 3 D) | 2 

= 4 (\D\ 2 - \B\ 2 ) + is 23 4 (afD - a*f*D* - D + D* 

- afB + a*f*B* + B - B*) + 44 (\a\ 2 - af - a* f* 
-BD* - B*D - \a\ 2 + af + a*f* + BD* + B*D) 
-is 3 23 c 23 (afB - a*f*B* - B + B* - afD + a*f*D* 
+D-D*)+st 3 {\B\ 2 -\D\ 2 ), 

AP^ ~ |S /i23 | 2 - |S /)32 | 2 = | - S23 c 23 (a* - /) - i/ (4^ " 4^) I 2 

- | - s 23 c 23 (a* -f)+ if {s 2 23 B - c\ 3 D) \ 2 

= 4 (|P| 2 - |P| 2 ) + is 23 4 (afB - a*f*B* - B + B* 

- afD + a*f*D* + D - D*) + 44 (\a\ 2 - af - a* f* 
-BD* - B*D - \a\ 2 + af + a*f* + BD* + B*D) 
-is 3 23 c 23 (afD - a*f*D* - D + D* - afB + a*f*B* 

+ B - B*) + 4 (|Pf - |P| 2 ) , 
AP r c r PT ~ |S /i33 | 2 - fe 3 | 2 = \s 2 23 a* + 4/ + is 23 c 23 / (P + D) | 2 

- I4«* + c 2 23 f + is 23 c 23 f (B + D) | 2 

= is 23 4 (B - B* + D- D* - B + B* - D + D*) 
+44 («/ + "7* + W + l^l 2 + BD* + B*D 
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- af-a*f* - \B\ 2 - \D\ 2 - BD* - B*D) 
+is 3 23 c 23 (afB - a*f*B* + afD - a*f*D* - afB 
+ a*f*B* - afD + a*f*D*) + s A 23 (|a| 2 - |a| 2 ) . (54) 

Note that AP e ^ PT is the only CPT probability difference that is uniquely determined by 
the (l,2)-subsector of the full three flavor neutrino evolution, see the explicit expressions 
of the evolution operators for neutrinos and antineutrinos [Eqs. flA35|) and (|A40j) ]. Thus, 
it is completely independent of the CP violation phase 8qp as well as the fundamental 
neutrino parameters Am^ ~ Am^, #13, and #23- 

Now, using conservation of probability, i.e., Eqs. (jHJ) - (fT3|) . we find the relations 

E ^P?I T = \C\>- \A\ 2 = 0, (55) 

a=e,/i,r 

E AP - T = \ A \ 2 ~ \ C \ 2 = °> ( 56 ) 

a>=e,{i,T 

E AP £ PT + E Ap ™ PT = E AP <T + E AP « CPT 

a=e,//,T a=e,/j,,r a=e,/i,,r a=e,fi,r 

= \B\ 2 + \D\ 2 - \B\ 2 - \D\ 2 = 0. (57) 
Thus, the CPT probability differences can be further simplified and we obtain 

APCPT „ |£|2 _ | /? |2 ) (5g) 

AP e c / T ^ 4 (|/3| 2 - |^| 2 ) - 2c 23S23 S (0/C7 - 0/^*) , (59) 

AP e c T PT ~ s 2 23 (\P\ 2 - \P\ 2 ) + 2c 23 s 23 % (13 fC - Pf*A*) , (60) 

Ai£ PT * 4 (\{3\ 2 - \f3\ 2 ) - 2c 23 s 23 % (Pf*A* - (3fC) , (61) 

AP T c e PT ~ s 2 23 (\(3\ 2 - \(3\ 2 ) + 2c 23 s 23 % ({3f*A* - PfC) , (62) 

where we have only displayed the CPT probability differences AP e ^ PT , AP e ^ PT , AP e ^ PT , 
AP^f T , and AP^, PT , since the remaining ones are too lengthy expressions and not so illu- 
minating. In the following, we will restrict our discussion only to those CPT probability 
differences displayed above. Furthermore, from the definition of the parameters a and b in 
Eq. (JA13J) . we can conclude that \b/a\ cx S 2 / A 2 = (Am^/Am^) 2 , and thus, the ratio \b/a\ 
is small, since Am 21 <C Am^. In Ref. [43], it has been shown that 



l«Mo)/WMo)l ~ \Ip,t (t,to)/I* a *, to (t,t )\ ~ ^ 2 /A 2 , 
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and therefore, it also holds that 



wmo)/wmo)i ~ m 



(Mo)/WMo)|~<* 2 /A 2 . 



Thus, the contributions of the integrals Ipj(t,to), I^ tQ {t,t ), Ip,t{b, to), and I% tQ (t,to) are 
suppressed by a factor of S 2 /A 2 in Eqs. f)59jl - Using this to reduce the arguments of 

the imaginary parts in Eqs. (JoTJj) - further, we obtain the following 



D. Examples of Matter Density Profiles 

We have now derived the general analytical expressions for the CPT violation probability 
differences. Next, we will calculate some of the CPT violation probability differences for 
some specific examples of matter density profiles. This will be done for constant matter 
density and step-function matter density. 

1. Constant Matter Density Profiles 

The simplest example of a matter density profile (except for vacuum) is the one of constant 
matter density or constant electron density. In this case, the matter potential is given by 
V(t) = V = const. Vi. Furthermore, if the distance between source and detector (i.e., the 
neutrino propagation path length or baseline length) is L and the neutrino energy is E u , 
then we can define the following useful quantities 



(3fC-(3f*A* ~ /3/a*I^ - 



(63) 
(64) 




(68) 



(65) 



(66) 



(69) 



(67) 
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(70) 



where 5 = ^,™ 21 . A 



2E V ' 



Am 31 
2E V 



Am: 



:!_> 



2£„ ' 



App.EJ) 



and 9 12 is the solar mixing angle. Then, we have (see 



(71) 
(72) 
(73) 
(74) 
(75) 
(76) 



a(t,0) 


= cos ojt + i cos 29 m 


sin tut, 


a(t,0) 


= cos tot + i cos 29 m 


sin cut, 


p(t,o) 


= — isin2# m sinu;t, 






= — i sin 2# m sin u)t, 




f(t,o) 


= e" iAt , 




f(t, 0) 


= e- iAt , 





where < t < L, which yield 



sin 2 29 m s 2 — sin 2 29 m s 

2 



-2 



s 2 c- S 2 ' S 

A, 12 c 12 w I 9 



UJ 2 



(77) 



3 (/3fC - Pf*A*) ~ s 12 c 12 s 13 5 (A - s 2 2 <5) | - ^\ cos(5cp 

(As — Q sin AL)s (As — cusin AL)s 



+ (A-c 2 2 5) 
+ (A-c 2 2 5) 



^2^2 _ A 2) 

(cos AL — c)s 



oj 2 {uj 2 — 
(cos AL — c)s 



A 2 ) 
sin Sqp 



cos 5, 



CP 



uo{uo 2 -K 2 ) lu(lu 2 -A 2 ) 
3 - PfC) ~ s 12 c 12 s 13 5 (A - s 2 2 5) | ^ - ^ cos5 C p 

(As — sin AL)s (As — u; sin AL)s 



(78) 



+ (A-c 2 2 5) 
- (A -c\ 2 5) 



u) 2 (u 2 -A 2 ) u 2 (u 2 -A 2 
(cosAL — c)s (cosAL — c)s 



cos 5, 



CP 



sin 5, 



CP 



(79) 



cu(cj 2 -A 2 ) cj(cj 2 -A 2 ) 

where s = sintuL, s = sinuL, c = coscjL, and c = cosu)L. Note that the only difference 
between the imaginary parts in Eqs. (fTSj) and (|T9~j) is the signs in front of the sin Sqp terms, 
i.e., applying the replacement Sep — > Sep, one comes from Q(j3fC — j3f*A*) to Q(/3f*A* — 
(3fC), and vice versa. Thus, inserting Eqs. (|77jl - (J79J) into Eqs. (J58j) - (JB2J), we obtain the 
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CPT probability differences in matter of constant density as 
ApCPT _ 2 2 r2 / sm 2 ^L sin 2 ^ 



UJ< 



UJ< 



A pCPT 



_ ,2 2 2 .2 
— i, 12 c 12 c 23 u 



sin 2 wL sin 2 tt;L\ 



x 



/ sin 2 ujL sin 2 u;L 



V ^ 2 



-^2 — J - 2S12C12S13S23C235 (A - s\ 2 5) 



uj' 



cos 5 



CP 



+ (A 
+ (A 



(As — u sin AL)s (As — c<j sin AL)s 



uj 2 (uj 2 - A 2 ) 



^2(^2 _ A 2) 



COS 5, 



CP 



(cos AL — c)s (cos AL — c)s 



AP 



CPT 



~ s 2 r 2 s 2 tf 2 

— 6 12 c 12 5 23° 



uj(uj 2 -A 2 ) uj(uj 2 -A 2 ) 
sin 2 cuL sin 2 CuL 



sin <5, 



CP 



x 



( sin 2 c<)L sin 2 ujL 



cos <5, 



+ 2S12C12S13S23C235 (A - s\ 2 5) 



CP 



+ (A 
+ (A 



(As — uj sin AL)s (As — uj sin AL)s 

w 2 (w 2 -A 2 ) cu 2 (cu 2 -A 2 ) 
(cos AL — c)s (cos AL — c)s 



cos <5, 



CP 



AP, 



CPT 

fie 



~ s 2 r 2 r 2 tf 2 

— i 12 c 12 c 23 w 



cu(cu 2 -A 2 ) cj(cj 2 -A 2 ) 
sin 2 u;L sin 2 ujL 



sin 5, 



CP 



X 



(sin 2 c<)L sin 2 cuL 
—2 



UJ 

cos 5qp 



2S12C12S13S23C235 (A - s 2 l2 S) 



+ (A 
- (A 



(As — uj sin AL)s (As — c<j sin AL)s 



cu 2 (cu 2 - A 2 ) 



u; 2 (u; 2 



cos 5, 



CP 



(cos AL — c)s (cos AL — c)s 



AP, 



CPT 



~ s 2 r 2 s 2 () 2 

— 6 12 c 12 6 23° 



uj(uj 2 -A 2 ) uj(uj 2 -A 2 ) 
sin 2 ujL sin 2 u)L 



sin 5, 



CP 



w 2 



+ 2S12C12S13S23C235 (A - s 2 2 5) 



x 



/ sin 2 ujL sin 2 c<jL 



uj' 



cos 5cp 



+ (A 
- (A 



(As — uj sin AL)s (As — c<j sin AL)s 



cj 2 (cu 2 - A 2 ) 
(cos AL — c)s 

wto 2 - A 2 ) 



cj 2 (cj 2 - A 2 ) 
(cos AL — c)s 



cos 5, 



CP 



cj(cj 2 - A 2 ) 



sin 5, 



CP 
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It is again interesting to observe that the CPT probability difference AP e ^ PT contains only a 
constant term in the mixing parameter Sqp, i.e., it is independent of the CP violation phase 
Sep, whereas the other CPT probability differences contain such terms, but in addition also 



sin5 CP and cos^cp terms (in the case of CP violation, see, e.g., Ref. |81(). Naively, one would 
not expect any sin<5cp terms in the CPT probability differences, since they do not arise in 
the general case of the T probability difference as an effect of the presence of matter, but 
are there because of the fundamental T violation that is caused by the CP violation phase 
5qp [43]. However, since constant matter density profiles are symmetric with respect to the 
baseline length L, the T violation probability difference is anyway actually equal to zero in 
these cases. Furthermore, we note that if one makes the replacement Sep — > — Sqp, then 
APC PT -> APC PT and ApC pT -> AP r c e PT and in the case that 5 CP = one has ApC pT = 
AP^f T and AP e ^ PT = AP^, PT . Moreover, in the case of degenerate neutrino masses mi = m 2 
or for extremely high neutrino energies, E u — > oo, the quantity S = A ^ 21 goes to zero and so 
do (3 and f3 (see the second point in the discussion at the end of App. |X] about the relation 
between Q and Q), which in turn means that the CPT probability differences in Eqs. (I58j) - 
(JfiSJl as well as in Eqs. flHDl) - flEU will vanish, i.e., AP^ PT -> when S -> 0. This can be 
understood as follows. In the case when Am|i <C Am^ (i.e., S <C A) or in the limit 5^0, 
we have that the neutrino mass hierarchy parameter £ = = ^ also goes to zero. If 



£ -> 0, then P ee -> 1 - sin 2 2g 13 sin gP AL , where C 13 = yW 29 13 + (f V - cos2# 13 ) 2 . Now, 
since we have only calculated the CPT probability differences to first order in perturbation 
theory in the small leptonic mixing angle 9 13 (see App. EJ, we have that P ee — > 1 when 
£ — > 0. Using P ee = 1 together with the unitarity conditions and (0), we find that 
Phij, = Ptt = 1 and P eAt = P^e = P^t = P T (i = 0, which means that neutrino oscillations 
will not occur in this limit. A similar argument applies for the case of antineutrinos. Thus, 
the CPT probability differences AP^ T — > up to first order in perturbation theory in 9 13 
when 5^0 (i.e., when S is completely negligible compared with A). Therefore, there are 
no extrinsic CPT violation effects up to first order in # 13 when 5^0. 

In the low-energy region V < S <C A, we find after some tedious calculations that 

nV ,rp r, / SL SL\ SLV 

APtf 1 ~ 8s 2 12 c 2 l2 cos 26»i2 [SL cos — - 2 sin — sin — — 

\ 2 2 J 2 

+ O ((V/Sf) , (85) 

. „PPT 9 9 9 I „ 5-^ 5L\ (5L V 

AP^ 1 ~ -8sf 2 q 2 Q3 cos 20i2 I 5Lcos — - 2 sin — ] sin— — 



2 2 ) 2 S 
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16si 2 cf 2 Si3S 2 3C23 cos 5qp cos 29 12 ( SL cos ^ 2 sin ] sin 



SL 

Y 



SL\ . SLV 
2 i ' 2~T 



16Si 2 Ci2Si 3 S23C23Sin5cP< COS 26»i2 



/. SL SL . 

x I oL cos — 2 sin — ) — sin SL 

+ o {{V/Sf) , 



SL cos SL — cos AL 



SL A I V 

+ oL sin — sin AL > — 

2 o 



AP° PT ~ -8si 2 Ci 2 S2 3 cos 26»i2 ( SL cos — - 2 sin 



5L 

y 



5L\ 5LV 

T sm TT 



+ 16si2cf 2 si3-S23C 2 3 cos 5cp cos 26*12 ( <5L cos — 2 sin — ] sin — — - 

2 2 J 2 o 



+ 16Sl2Cl 2 Sl3S23C23Sin5 C PS COS 26*12 



SL cos SL — cos AL 



cr . SL . . T \V 
+ oL sin — sin AL > — 
2 o 



x oL cos ism — — sin oL 

V 2 2) 

+ ^ ((VV5) 3 ) , 

. „ rPT ^999 /c-^ SL „ . <5L\ . 5L V 

AP^/ 1 ~ -8sf 2 ci 2 C23 cos 26»i2 ULcos— -2 sin — 1 sin— — 



— 16si2cf 2 -Si3S23C23 cos Sep cos 26*12 ( SL cos — 2 sin — ] sin — — 



SL 

Y 



2 5 
SL 
2 



S_LV_ 

YJ 



+ 16Sl 2 Ci 2 Sl3S23C23Sin(5cP< COS 26*12 



x ( SL cos ^ — 2 sin ^ ] - sin SL 



SL cos 5L — cos AL 



cr . SL . AT \V 
+ oL sin — sin AL > — 
2 | 5 



AP T c e PT ~ -8si 2 Ci 2 S2 3 cos 26*12 ( SL cos ^- - 2 sin 



<5L 

y 



5L\ <5LV 

T sm TT 



+ 16si 2 Ci 2 Si 3 s 2 3C23 cos S C p cos 26*i 2 I SL cos — 2 sin ^ J sin 

2 2 I 2 o 



16si 2 Ci 2 Si3S 2 3C9 3 sin£ C p<( cos26*i 2 



x j 5Lcos^ — 2sin^ ) - sin e)! 



SL cos SL — cos AL 



cr . SL . . T \ V 
+ oL sin — sin AL > — 
2 | S 



+ O ({V/ 8?) . 



(86) 



(87) 



(88) 



(89) 



Note that there are, of course, no terms in the CPT probability differences that are constant 
in the matter potential V, since in the limit V — > 0, i.e., in vacuum, the CPT probability 
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differences must vanish, because in vacuum they are equal to zero [cf., Eq. flU]) ]. Further- 
more, we observe that the leading order terms in the CPT probability differences are linear 
in the matter potential V, whereas the next-to- leading order terms are cubic, i.e., there are 
no second order terms. However, we do not show the explicit forms of the cubic terms, since 
they are quite lengthy Actually, for symmetric matter density profiles it holds that the 
oscillation transition probabilities in matter for neutrinos and antineutrinos, P(u a — > vp; V) 
and P{vn — > up;V), respectively, are related by P{y a — > vp;V) = P{pp — > u a ; — V) [see 
Ref. [60j and Eqs. pUj) and (jHJ)]. Hence, in this case, the CPT probability differences 
AP a C J T (V) = ?K -> up; V) - P{yp -> u a , V) = P(u a -> up; V) - P(u a -> up; -V) = 
f(V) — /(— V) are always odd functions with respect to the (symmetric) matter potential 
V, since AP^(-V) = /HO ~ f(V) = - 



TtP 

Introducing the Jarlskog invariant 



asllaai 



f(V) - f(-V)] = -AP^(V) [9 



J = S12C12S13C13S23C23 sin 5 C p ^ S12C12S13S23C23 sin 5 C p, (90) 
we can, e.g., write the CPT probability difference AP e ^ PT as 



AP£ PP * -c4AP e ° e PT 



16c^ 2 cos 2^12 J cot <5cp ( 5i cos — 2 sin — ) sin — — 



- 16 J< cos 2d 



12 



5L\ 5LV 

2 sin — sin - 

2 2 7 2 5 

, , , . , 5L SL 

oL cos oL — cos AL oL cos 2 sin — 

2 2 



+ 5Lsin^sinAL + ((V/5) 3 ) . (91) 



— sin5L 

In the case of maximal solar mixing, i.e., if the solar mixing angle #12 = \, then we have 

AP e c e PT ~ 0, (92) 
which is also obtained using Eq. (|50jl . and 



AP e c / T ~ -16 J5L sin ^ sin ALy ~ -AP£ PT , (93) 



where in this case J = |si 3 S23C 2 3 sin Sep- Thus, we would not be able to observe any extrinsic 
CPT violation in the u e — > u e and u e —> u e channels. However, it would still be possible to 
do so in the u e — > u^ and v e — > channels. Furthermore, note that if in addition <5cp = 0, 
then also APj?, PT and AP£, PT vanish, since J oc sin^cp. 



18 



2. Step-function matter density profiles 



Next, we consider step-function matter density profiles, i.e., matter density profiles con- 
sisting of two different layers of constant densities. Let the widths of the two layers be L\ 
and L 2 , respectively, and the corresponding matter potential V\ and V 2 . Furthermore, we 
again let E v denote the neutrino energy. Similar to the constant matter density profile case, 
we define the quantities 



OJi = 

COi = 

Aj = 

Ai = 

u m,i — 

ft = 



Oft Vi 
cos 20i2 — r 

o 



+ sin 2 20 



12) 



- \ I cos 20i2 + y 



+ sin 2 20 



12- 



1 /5cos20i 2 -^ 

- arccos 

2 V 

f /5cos20i2 + V- 

- arccos 

2 V 2cui 



with i = 1,2 denoting the two different layers, where again 5 



Am'i 



A 



Ami 



(94) 

(95) 
(96) 
(97) 
(98) 
(99) 
and 0i2 is 



2E V ' " — 2E y ' 

the solar mixing angle. We divide the time interval of the neutrino evolution into two parts: 
< t < L\ and L\ < t < L, where L = L\ + L 2 . In the first interval, the parameters a, 
a, (3, (3, f, and / are given by the well-known evolution in constant matter density, i.e., by 



Eqs. ffTTj) - (fTtjj) with the replacements u> — > uj\, uj — > u>i, A — > A 1} A — > Aj, m 
#m — ► ®m i) whereas in the second interval, they are given by 



m i, and 



a (Mo) 
a(t,t ) 

Pit, t ) 
P(t,t ) 
f(t,t ) 
/(Mo) 



= ci4 - sis 2 cos(20 mi i - 20 mi2 ) 

+ i(sic' 2 cos20 mj i + s' 2 Ci cos29 m 2 ), 

= c~ic 2 - si s 2 cos(20 mi i - 20 m _ 2 ) 

+ i(s"ic 2 cos 20 m i + s' 2 c\ cos 26 m 2 ), 

= sis' 2 sin(20 m) i - 20 mi2 ) - i(sic 2 sin20 mi i + s 2 c x sin 20 m)2 ), 

= sis 2 sin(20 mi i - 20 mi2 ) - i(sic 2 sin ^m,i + s' 2 c~i sin 29 m>2 ), 

_ e -i[Ai(L 1 -t )+A a (i-i 1 )] 

_ -i[Ai(Li-to)+A 2 (t-ii)l 



(100) 

(101) 
(102) 
(103) 
(104) 
(105) 
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where Sj = smc^Li, q = coscUjLj, Sj = sinc^iLj, and q = coscjjLj with i = 1,2 and 
s 2 = sinter, c' 2 = coso^t, s 2 = sin^r, and c 2 = cosa) 2 r with r = t — L\. 

Now, we will take a quick look at the general way of deriving expressions for the CPT 
probability differences for the step-function matter density profile. However, in this case, 
the derivations are quite cumbersome and we will only present the results for the CPT 
probability difference AP e ^ PT . 

Similar to the case of constant matter density, we obtain the CPT probability difference 
AP^ for step-function matter density profiles as 

AP° P T ~ (sic 2 cos 20 m ,i + s 2 c l cos29 ni)2 ) 2 + [c\c 2 - Sis 2 cos2(0 m) i - 9 m)2 )} 2 
- (sic 2 cos 29 mA + s 2 ci cos 2# mi2 ) 2 + [c x c 2 - Sxs 2 cos2(9 m>1 - m)2 )] 2 . 

(106) 



In the low-energy region V\ t2 < 5 <C A, we find that 
AP£ PT ~ 8s 2 c 2 o cos20 



V x V 2 \ 5(L 1 + L 2 ) 
6 L x — + L 2 — cos 



6 6 J 2 
n fVx . 5L 1 5L 2 V 2 . 5L 2 5 LA] . 5{Lx + L 2 ) 

~ 2 [j sm — cos — + T sm T cos T ) \ sm 2 

+ O ((Fx/5) 2 , (V 2 /5) 2 , VxV 2 /5 2 ) . (107) 

One observes that the CPT probability difference AP e ^ PT is completely symmetric with 
respect to the exchange of layers 1 and 2. Furthermore, in the limit Vx, 2 — > V and Lx j2 — > L/2, 
one recovers the CPT probability difference for constant matter density (as one should), see 
Eq. (JEU). 

III. NUMERICAL ANALYSIS AND IMPLICATIONS FOR NEUTRINO OSCIL- 
LATION EXPERIMENTS 

In general, the three flavor neutrino oscillation transition probabilities in matter P a p = 
P(v a — > vp) are complicated (mostly trigonometric) functions depending on nine parameters 

P a /3 = P*/3(Am 21 , Anij!, 9 12 , 9 13 , 9 23 , 5 C p; E v , L, V(L)), a, (3 = e, fi,r, (108) 

where Am 21 and Am^ are the neutrino mass squared differences, 6 I 12 , 9x3, 9 23 , and 5qp are 
the leptonic mixing parameters, E v is the neutrino energy, L is the baseline length, and 
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TABLE I: Present values of the fundamental neutrino parameters. 



Parameter Best-fit value Range References 





7.1 • 10" 5 eV 2 


~ (6 4- 9) • io- 5 eV 2 (99.73 % C.L.) [27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 42, ,100, 101} 


\Am 2 31 \ 


2.5 • 10" 3 eV 2 


(1.6 4- 3.9) ■ IO" 3 eV 2 (90 % C.L.) [23, 24, .25, .26] 


#12 


34° 


27° + 44° (99.73 % C.L.) [27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 42, ,100, 101} 


#13 




9.2° (90 % C.L.) [102, .103, .104] 


#23 


45° 


37° -=- 45° (90 % C.L.) [23., 24, 25., 26J 


<5cp 




[0,2vr) 



finally, V(L) is the matter potential, which generally depends on L. Naturally, the CPT 
probability differences depend on the same parameters as the neutrino oscillation transition 
probabilities. The neutrino mass squared differences and the leptonic mixing parameters are 
fundamental parameters given by Nature, and thus, do not vary in any experimental setup, 
whereas the neutrino energy, the baseline length, and the matter potential depend on the 
specific experiment that is studied. 

The present values of the fundamental neutrino parameters are given in Table |U These 
values are motivated by recent global fits to different kinds of neutrino oscillation data. 
All results within this study are, unless otherwise stated, calculated for the best-fit values 
given in Table Furthermore, we assume a normal neutrino mass hierarchy spectrum, i.e., 
Amjj <C A m 3i with Amjj = +2.5 ■ 10~ 3 eV 2 . For the leptonic mixing angle #13, we only 
allow values below the CHOOZ upper bound, i.e., sin 2 29i 3 < 0.1 or #13 < 9.2°. For the 
CP violation phase, we use different values between and 2tt, i.e., 5qp € [0, 27r). Note that 
there is no CP violation if 5qp £ {0, rr}, whereas the effects of CP violation are maximal if 
<5cpG{|,f}. 

As realistic examples, let us now investigate the effects of extrinsic CPT violation on the 
transition probabilities for neutrino oscillations in matter for various experiment in which 
the neutrinos traverse the Earth. Such experiments are e.g. so-called long baseline exper- 
iments, atmospheric and solar neutrino oscillation experiments. In some analyses of these 
experiments, the Preliminary Reference Earth Model (PREM) matter density profile jl05| 
has been used, which has been obtained from geophysics using seismic wave measurements. 
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However, the (mantle-core-mantle) step-function m atter density profile 135] is an excellent 



approximation to the PREM matter density profile 106], whereas the constant matter den- 
sity profile serves as a very good approximation to long baseline experiments that have 
baselines that do not enter the core of the Earth. Thus, we use these approximations for 
our calculations. 

The equatorial radius of the Earth and the radius of the core of the Earth are ~ 
6371 km and r ~ 3486 km, respectively, which means that the thickness of the mantle of the 
Earth is Rq — r ~ 2885 km. From the geometry of the Earth, one finds that the relation 
between the maximal depth of the baseline t and the baseline length L is given by 



L 2 



£ = R (B -^Rl-— (orL = 2 v /£(2 J R e -£)). (109) 

Hence, in order for the neutrinos also to traverse the core of the Earth, i.e., I > R§ — r, the 
baseline length needs to be L > 10670 km. This means that for experiments with baseline 
lengths shorter than 10670 km, we can safely use the constant matter density profile. For 
"shorter" long baseline experiments (L < 3000 km) we use the average matter density of the 
continental Earth crust, p crus t — 3g/cm 3 , whereas for "longer" long baseline experiments 
(3000 km < L < 10670 km) we use the average matter density of the mantle of the Earth, 
Pmantie — 4.5g/cm 3 . Furthermore, the matter potential V = V(L) expressed in terms of the 
matter density p = p(L) is given by 

V ~ -^G F —p ~ 3.78 ■ HT 14 eV • p[g/cm 3 ], (110) 
V 2 m N 

where p[g/cm 3 ] is the matter density given in units of g/cm 3 . 

Let us now investigate when it is possible to use the low-energy approximations for the 
CPT probability differences derived in the previous section. In these approximations, we 
have assumed that the matter potential V is smaller than the parameter S, i.e., V < 5 A. 
Now, the parameter 5 is a function of the neutrino energy E u : 



2 

3.55- 10- 5 eV ■ E^eV}- 1 , (111) 



Am 2i _ 



2E U 

where -E^feV] is the neutrino energy in eV. Thus, combining Eqs. (jllOj) and (jlll|) . we find 
that 

E u < 0.94 • 10 9 eV • pfg/cm 3 ] -1 , (112) 
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which means that for the continental Earth crust (p C rust — 3g/cm 3 ) the neutrino energy 
E v must be smaller than about 0.31 GeV in order for the low-energy approximations to be 
valid. 

In Table HTl we list several past, present, and future long baseline experiments of accel- 
erator and reactor types including their specific parameter sets for which we are going to 
estimate the extrinsic CPT violation effects. From the values of the neutrino energies given 
in this table we can conclude that the low-energy approximations for the CPT probabil- 
ity differences are applicable for the reactor experiments including the LSND accelerator 
experiment, but not for the accelerator experiments in general. 

Using the values of the fundamental neutrino parameters given in Table as well as the 
approximate values of the neutrino energy and baseline length for the different long baseline 
experiments given in Table HTl we obtain estimates of the CPT probability differences, which 
are presented in Table IIHI From the values in Table IIHI we observe that there are three 
different experiments with fairly large estimates of the CPT probability differences. These 
experiments are the KamLAND, BNL NWG, and NuMI experiments, which will later in this 
paper be studied in more detail. In general, there is a rather large discrepancy among the 
values coming from the numerical, analytical, and low-energy approximation calculations. 
This is mainly due to the oscillatory behavior of the CPT probability differences. Therefore, 
these values can change drastically with a small modification of the input parameter values. 
Thus, this can explain the somewhat different values of the different calculations. However, 
in most of the cases, the order of magnitude of the different calculations are in agreement. 
Note that for all reactor experiments the analytical and low-energy approximation estimates 
agree completely, since the neutrino energies are low enough for these experiments in order 
for the low-energy approximations to be valid. Moreover, we have calculated the CPT 
probability difference APf^" 1 for two potential neutrino factory setups using the analytical 
formula (}83|) . In general, these setups are very long baseline experiments that even penetrate 
the Earth's mantle in addition to the Earth's crust. For our calculations we used a constant 
matter density profile with p = p man tic — 4.5g/cm 3 . Furthermore, we chose the neutrino 
energy to be 50 GeV as well as the baseline lengths 3000 km and 7000 km, respectively. For 
these parameter values, we obtained AP^/ 1, ~ 3.0 • 10~ 5 (3000 km) and AP^ e PT ~ 1.8 • 10~ 5 
(7000 km). Thus, the extrinsic CPT violation is practically negligible for a future neutrino 
factory. 
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TABLE II: Accelerator and reactor long baseline experiments including measurable neutrino os- 
cillation channels, average neutrino energies (E u ), approximate baseline lengths (L) as well as 
references to the respective experiments. The CHOOZ, KamLAND, and Palo Verde experiments 
are reactor experiments, whereas the other experiments are accelerator experiments. Furthermore, 
the BooNE, MiniBooNE, CHOOZ, LSND, NuTeV, and Palo Verde experiments are sometimes 
called short baseline experiments. However, we will use the term long baseline experiments for all 
experiments in this table. 



Experiment 


Channels 


E v 


L 


References 


BNL NWG 


v u — 
/* 


> v e 


lGeV 


400 km, 2540 km 


ri07. 108. 1091 


1 

BooNE 4 


1 P 
l % ~ 


-> v e 

-> V e 


(0.5 "t" 1.5) GeV 


1 km 


[no] 


MiniBooNE j 




->■ V e 
->■ P e 


(0.5 -=- 1.5) GeV 


500 m 


PJUJ 


CHOOZ 




> 


~ 3MeV 


1030 m 


ri02. 103. 1041 


ICARUS | 


I ^ - 

I Vp, - 


->■ zy e 

-> V T 


17 GeV 


743 km 


[112, U3, 114J 


JHF-Kamioka j 


K P 


-> 

P 


(0.4 ~r~ 1.0) GeV 


295 km 


[115] 


I 

K2K <^ 

I 


I w P 

U,i - 


— > ZA, 
A 4 


1.3 GeV 


250 km 


[116,112] 


KamLAND 


v e - 


> 


~ 3MeV 


~ 180 km 


[411 


LSND | 


I v n - 

{ Vp - 


-> Z/ e 
-> P e 


48 MeV 


30 m 


[20, 21, 22] 


MINOS | 


I ^ - 

[ - 


-> Z/ e 
"> ^ 


(3 ~r~ 18) GeV 


735 km 


ri!8. 119. 1201 


NuMI I /ll j 


I ^ - 


-> f e 


1.4 GeV / 0.7 GeV 712km / 987km 


[121] 




t - 


-> P e 








NuTeV | 


[ Vp - 

{ ^ - 


-> Z/ e 

-> 9 e 


75 GeV, 200 GeV 


(915 -=- 1235) m 


[122] 


OPERA 


Vp - 


> i/ T 


17 GeV 


743 km 


[114, .123] 


Palo Verde 


V e - 


> 


~ 3 MeV 


750 m, 890 m 


[124 , .125 , .126 , .127] 
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TABLE III: Estimates of the CPT probability differences for the different long baseline experiments 
listed in Table ILTI The fundamental neutrino parameters used are: Arn^i = 7.1 • 10 -5 eV 2 , Am^ = 
2.5 • 1CT 3 eV 2 , 9i2 = 34°, 9 13 = 9.2°, (9 23 = 45°, and 5 C p = 0. Furthermore, we have used constant 
matter density profiles with p = 3g/cm 3 as approximations of the continental Earth crust. 
Experiment CPT probability differences 



Quantities Numerical Analytical Analytical (low-energy) 



BNL NWG 




0.010 


3.6- 10~ 4 


1.7- 10~ 6 


1 Jl \ l_j 1 > V V VI 


ApCPT 


03? 


1.2 • 10~ 3 


9 7 • 10 -3 


RnnNE 


ApCPT 


6 6-1 o -13 


5 1 • 10~ 14 


9 0-1 -17 


MiniRooNF 

IVlllllJJWvl > 1 J 


A pCPT 


4.1 • 10~ 14 


3.2 • 10~ 15 


—2.0 • 10~ 17 


CHOOZ 


APCPT 


-3.6 • 10~ 5 


-3.7 • 10~ 9 


-3.7- 10~ 9 


TCAPTTS 


ApCPT 
ApCPT 


4 0-1 -5 

— 3 R • 1 0~ 5 


3 1 • in -6 


4 1-1 -9 

L ±. ± J-\J 


XYdlIJ.lVJ.LVCl 


ApCPT 

ApCPT 
W 


3 8-1 -3 
—13-1 0~ 4 


9 9 • in -4 


5 0-1 o -7 


K9K 


ApCPT 

ApCPT 
W 


10-1 o -3 
—53-1 0~ 5 


7 9 • in -5 


19-1 -7 


KamLAND 


ApCPT 


—0 033 


—0 040 


—0 040 


LSND 


A pCPT 

fie 


4.8 • 10~ 15 


3.7 • 10~ 16 


1.9 • 10~ 18 


MINOS 


AP£ PT 


1.9 • 10" 4 
-1.1 • 10~ 5 


1.4- 10~ 5 


1.9 • 10~ 8 


NuMI I 


AP, CPT 


0.026 


-2.7 • 10~ 5 


6.2 • 10~ 6 


NuMI II 


AP, CPT 


2.6 • 10~ 3 


-2.4 • 10~ 4 


1.8 • 10" 4 


NuTeV 


AP£ PT 


1.6 • 10~ 18 


1.2 • 10~ 19 


-2.6 • 10~ 15 


NuTeV 


AP£ PT 


8.2 • 10" 20 


6.4- 10~ 21 


-1.5 • 10~ 15 


OPERA 


AP, CPT 


-3.8 • 10~ 5 






Palo Verde 


APCPT 


-1.2 • 10~ 5 


-1.1 • 10~ 9 


-1.1 • 10~ 9 


Palo Verde 


APCPT 


-2.2 • 10~ 5 


-2.1 • 10~ 9 


-2.1 • 10~ 9 



Next, in Fig. ^ we plot the CPT probability differences AP e ^ PT and AP^ PT as functions 
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of the neutrino energy E u for three different characteristic baseline lengths: 1 km, 250 km, 
and 750 km. From these plots we observe that the CPT probability differences increase with 




FIG. 1: The CPT probability differences APCPT 

and AP^ e PT plotted as functions of the neutrino energy E v . The baseline 
lengths used are: 1 km (left column), 250 km (middle column), and 750 km (right column) with p = 3g/cm 3 . Do t ted c urve s 



correspond to numerical calculations using the evolution operator method and Cayley-Hamilton formalism 
whereas solid curves correspond to analytical calculations using Eqs. 1801 and 1831 . The fundamental neutrino parameters used 
are: C\m\ x = 7.1 ■ 10" 5 eV 2 , Am^ = 2.5 ■ 10" 3 eV 2 , 6»i 2 = 34°, 9 13 = 9.2°, 6» 2 3 = 45°, and <5 C p = 0. 



increasing baseline length L. Furthermore, we note that for increasing neutrino energy E u 
the extrinsic CPT violation effects disappear, since the CPT probability differences go to 
zero in the limit when E v — > oo. We also note that AP e ^ PT and AP^, PT are basically of 
the same order of magnitude. In this figure, the numerical curves consist of a modulation 
of two oscillations: one slow oscillation with larger amplitude and lower frequency and 
another fast oscillation with smaller amplitude and higher frequency. On the other hand, the 
analytical curves consist of one oscillation only and they are therefore not able to reproduce 
the oscillations with smaller amplitudes and higher frequencies. However, the agreement 
between the two curves are very good considering the oscillations with larger amplitudes and 
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lower frequencies. In principle, the analytical curves are running averages of the numerical 
ones, and in fact, the fast oscillations cannot be resolved by any realistic detector due to 
limited energy resolution making the analytical calculations excellent approximations of the 
numerical ones. 

Let us now investigate some of the most interesting experiments in more detail for which 
the extrinsic CPT violation effects may be sizeable. In Fig. we plot the CPT probability 
difference AP e ^ PT as functions of both the neutrino energy E v and the baseline length L 
centered around values of these parameters characteristic for the KamLAND experiment. 
We observe that for neutrino energies around the average neutrino energy of the KamLAND 



0.04 




,6 1Q 7 1Q 8 100 150 200 250 300 

E v [eV] L [km] 



FIG. 2: The CPT probability difference AP e ^ PT for the KamLAND experiment. The left-hand side plot shows its dependence 
on the neutrino energy E„, whereas the right-hand side plot shows its dependence on the baseline length L. The solid and 
dotted curves are analytical and numerical results, respectively. The diamonds ('©') indicate the central values of the KamLAND 
experiment. The parameters used are the same as for Fig. Ill 

experiment the CPT probability difference AP e ^ PT could be as large as 3 % - 5 % making 
the extrinsic CPT violation non- negligible. This means that the transition probabilities 
P( v e ~^ v e) and P{p e —>■ u e ) are not equal to each other for energies and baseline lengths 
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typical for the KamLAND experiment. Thus, if one would be able to find a source of electron 
neutrinos with the same neutrino energy as the reactor electron antineutrinos coming from 
the KamLAND experiment, then one would be able to measure such effects. Furthermore, 
for the KamLAND experiment the agreement between the analytical formula (J8Uj) and the 
low-energy approximation (J85)) is excellent, i.e., it is not possible to distinguish the results 
of these formulas from each other in the plots. 

Next, in Figs. H3 - El we present some plots for the topical accelerator long baseline 
experiments BNL NWG, JHF-Kamioka, K2K, and NuMI, which have approximately the 
same neutrino energies, but different baseline lengths. In these figures, we plot the CPT 
probability difference AP^, PT as functions of the neutrino energy E u and the baseline length 
L as well as the neutrino energy E v for three different values of the CP violation phase 
<5cp corresponding to no CP violation (5cp = 0), "intermediate" CP violation (5cp = f ), 
and maximal CP violation (<5cp = §), respectively. We note that in all cases the low- 
energy approximation curves are upper envelopes to the analytical curves. Furthermore, 
we note that the CPT probability difference AP^ T is larger for long baseline experiments 
with longer baseline lengths and it does not change radically for different values of the CP 
violation phase Sep- 




FIG. 3: The CPT probability difference AP^PT for the BNL nw q experiment (baseline length: 2540 km). The left-hand 
side plot shows its dependence on the neutrino energy E v (solid curve = analytical calculation, dashed curve = low-energy 
approximation), the middle plot shows its dependence on the baseline length L (solid curve = analytical calculation, dashed 
curve = low-energy approximation), and the right-hand side plot shows the dependence on E v for three different values of <5cp: 
(solid curve), 21 (dotted curve), and ^ (dashed curve). The diamonds ('o') indicate the central values of the BNL NWG 
experiment. The other parameters used are the same as for Fig. HI 
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FIG. 4: The CPT probability difference AP / J e PT for the JHF-Kamioka experiment. The left-hand side plot shows its 
dependence on the neutrino energy E v (dotted curve = numerical calculation, solid curve = analytical calculation, and dashed 
curve = low-energy approximation), the middle plot shows its dependence on the baseline length L (solid curve = analytical 
calculation), and the right-hand side plot shows the dependence on E v for three different values of 5cp : (solid curve), ^ 
(dotted curve), and J (dashed curve). The diamonds ('o') indicate the central values of the JHF-Kamioka experiment. The 
other parameters used are the same as for Fig. Q 




FIG. 5: The CPT probability difference AP° PT for the K2K experiment. The left-hand side plot shows its dependence on 
the neutrino energy E y (dotted curve = numerical calculation, solid curve = analytical calculation, and dashed curve = low- 
energy approximation), the middle plot shows its dependence on the baseline length L (solid curve = analytical calculation), 
and the right-hand side plot shows the dependence on E v for three different values of <5cp : (solid curve) , (dotted curve) , 
and 2- (dashed curve). The diamonds ('o') indicate the central values of the K2K experiment. The other parameters used are 
the same as for Fig. HI 

Finally, in Fig. we present numerical calculations shown as density plots of the CPT 
probability differences for neutrinos traversing the Earth, which are functions of the nadir 
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FIG. 6: The CPT probability difference AP° e PT for the NuMI Phase II experiment. The left-hand side plot shows its 
dependence on the neutrino energy E v (dotted curve = numerical calculation, solid curve = analytical calculation, and dashed 
curve = low-energy approximation), the middle plot shows its dependence on the baseline length L (dotted curve = numerical 
calculation, solid curve = analytical calculation, and dashed curve = low-energy approximation), and the right-hand side plot 
shows the dependence on E v for three different values of <Scp : (solid curve), ^ (dotted curve), and ^ (dashed curve). The 
diamonds ('©') indicate the central values of the NuMI Phase II experiment. The other parameters used are the same as for 
Fig.Hl 

angle h and the neutrino energy E v . The numerical calculations are based on the evol ution 



operator method and Cayley-Hamilton formalism introduced and developed in Refs. 



129 



128, 



1301 ] and the parameter values used are given in the figure caption. The nadir angle 
h is related to the baseline length L as follows. A nadir angle of h = corresponds to 
a baseline length of L = 2/2®, whereas h = 90° corresponds to L = 0. As h varies from 
to 90°, the baseline length L becomes shorter and shorter. At an angle larger than 
hn = arcsin ~ 33.17°, the baseline no longer traverse the core of the Earth. The CPT 
probability differences in Fig. [7| might be of special interest for atmospheric (and to some 
extent solar) neutrino oscillation studies, since the plots cover all nadir angle values and 
neutrino energies between 100 MeV and 100 GeV. We note from these plots that for some 
specific values of the nadir angle and the neutrino energy the CPT violation effects are rather 
sizeable. 
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FIG. 7: The CPT probability differences AP^p T (a, f3 = e,fi, t) plotted as functions of the nadir angle h and the neutrino 
energy E v . The different plots show: AP£ PT (upper-left), AP^ T (upper- middle), AP^ PT (upper-right), AP^ e PT (middle- 
left), APCPT (middle- middle), AP^PT ( m iddle-right), AP T c e PT (down- left), ApC pT (down-middle), and ApC. PT (down-right). 
The fundamental neutrino parameters used are: Am^ = 7.1 • 10 _5 eV 2 , Am^ = 2.5 ■ 10 _3 eV 2 , 612 = 34°, #13 = 9.2°, 
$23 = 45° , and <5cp = 0. Furthermore, we have used the mantle-core-mantle step- function approximation of the Earth matter 
density profile. 

IV. SUMMARY AND CONCLUSIONS 

In conclusion, we have studied extrinsic CPT violation in three flavor neutrino oscil- 
lations, i.e., CPT violation induced purely by matter in an intrinsically CPT-conserving 
context. This has been done by studying the CPT probability differences for arbitrary mat- 
ter density profiles in general and for constant matter density profiles and to some extent 
step-function matter density profiles in particular. We have used an analytical approxima- 
tion based on first order perturbation theory and a low-energy approximation derived from 
this approximation as well as numerical calculations using the evolution operator method 
and Cayley-Hamilton formalism. The different methods have then been applied to a num- 
ber of accelerator and reactor long baseline experiments as well as possible future neutrino 
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factory setups. In addition, their validity and usefulness have been discussed. Furthermore, 
atmospheric and solar neutrinos have been studied numerically using a step-function matter 
density profile approximation to the PREM matter density profile. Our results show that 
the extrinsic CPT probability differences can be as large as 5 % for certain experiments, but 
be completely negligible for other experiments. Moreover, we have found that in general 
the CPT probability differences increase with increasing baseline length and decrease with 
increasing neutrino energy. All this implies that extrinsic CPT violation may affect neu- 
trino oscillation experiments in a significant way. Therefore, we propose to the experimental 
collaborations to investigate the effects of extrinsic CPT violation in their respective exper- 
imental setups. However, it seems that for most neutrino oscillation experiments extrinsic 
CPT violation effects can safely be ignored. 

Finally, we want to mention that in this paper, we have assumed that the CPT invariance 
theorem holds, which means that there will be no room for intrinsic CPT violation effects 
in our study, and therefore, the CPT probability differences will only contain extrinsic CPT 
violation effects due to matter effects. However, it has been suggested in the literature that 
there might be small intrinsic CPT violation effects in neutrino oscillations , which 
might be entangled with the extrinsic CPT violation effects. The question if such intrinsic 
and the extrinsic CPT violation effects could be disentangled from each other in, for example, 
realistic long-baseline neutrino oscillation experiments is still open J2] and it was not the 
purpose of the present study. Actually, this deserves an own complete systematic study. 
However, such a study would be highly model dependent, since intrinsic CPT violation 
is not present in the SM. Furthermore, it should be noted that in the above mentioned 
references, Refs. , the intrinsic CPT violation effects were only studied in neutrino 
oscillations with two flavors and not with three. 
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APPENDIX A: EVOLUTION OPERATORS 



Neutrino oscillations are governed by the Schrodinger equation [see Eq. 



i^K*)> = ^(*)K*)>- 



(Al) 



Inserting \u{t)) = S(t,to)\v(to)) [Eq. (|3*3j) ] yields the Schrodinger equation for the evolution 
operator 

(A2) 



i-S(Mo) = Jt(t)S(t,t ), 



which we write in flavor basis as 



i-s f (t,t ) = jr f (t)Sf(t,to). 



(A3) 



In what follows, we will assume that the number of neutrino flavors is equal to three, i.e., 
n = 3. Thus, the total Hamiltonian in flavor basis for neutrinos is given by 



je f (t) = H f + V f (t) = UH m U ] + V f (t), 



(A4) 



where 



H„ 




5 
A 



( V(t) 0^ 



and V f (t) 



\ 








are the free Hamiltonian in mass b asis a nd the matter potential in flavor basis, respectively, 



and U is the leptonic mixing matrix 



136j . Here 5 = 



2E V ' 



A = ^gi, and V(t) = V2G F N e (t) 
is the charged-current contribution of electron neutrinos to the matter potential, where 
Gf — 1.16639 • 10~ 23 eV -2 is the Fermi weak coupling constant and N e (t) = ^-p(t) is 
the electron number density with Y e being average number of electrons per nucleon (in the 
m N ~ 939.565330 MeV the nucleon mass, and p = pit) the matter density. 



Earth: Y P 



2)1 



The sign of the matter potential depends on the presence of neutrinos or antineutrinos. In 
the case of antineutrinos, one has to change the sign by the replacement V(t) — > —V(t). 
Thus, the the total Hamiltonian in flavor basis for antineutrinos is given by 



J?,(t) =H f - V,(t) = UH m U^ - V f (t). 



(A5) 



Decomposing U = O23U13O12 = O23U', we can write the total Hamiltonian in flavor basis 
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as 



J?f(t) = 23 



A) 0^ 



u' 



5 
A 



( V(t) (A 



\ 








T 23 = 2 zH{t)0^. 



(A6) 



Here we use the following parameterization for the orthogonal matrices O23 and 0\ 2 and the 
unitary matrix U\ 3 



(\ ^ 



Ok 



23 



C23 S23 
-S23 c 23 



( c 13 s 13 e- i5 cp\ / ■ 



17- 



13 



-si 3 e 



1 

i<5cP Q 



Cl3 



On 



12 



/ 



V 



Cl2 Si 2 

-S12 C12 
1 



where s a b = sm9 a b and c a b = cos^. Here #12, 6>i 3 , and # 23 are the ordinary vacuum mixing 
angles and 5 C p is the CP violation phase. This means that U is given by the standard 
parameterization of the leptonic mixing matrix and U' is given by 



U' = 



( 



V 



Cl 3 Ci2 Ci 3 Si2 si 3 e" 

-S12 C12 



-si 3 ci 2 e i<5cp s 13 s 12 e i5 ^ c 13 



-i5cp 



(A7) 



Inserting J0f(t) = 23 H(t)0 23 into the Schrodinger equation, we obtain 



i—S(t,t ) = H(t)S(t,t ), 



(A8) 



(A9) 



where S(t,t ) = 23 Sf(t, t )O 23 . Thus, the Hamiltonian H(t) can be written as 

f (? 13 sj 2 5 + sj 3 A + V(t) c 13 c 12 s 12 5 c 13 s 13 (A-sj 2 5)e-' lScp \ 
H(t) = ci 3 ci 2 si 2 5 cj 2 S s 13 c 12 s 12 e~ lScp 5 

yC 13 s n (A - s\ 2 8) e i<5cp -s 13 c 12 s 12 e iScp S s\ 3 s\ 2 5 + c 2 13 A 

Series expansions of s 13 and c i3 when 13 is small, i.e., s 13 = 9i 3 + 0(9f 3 ) and c 13 = l + 0(9f 3 ), 
gives up to second order in 13 



/ 



5s\ 2 + V(t) 



5c 12 s 12 



9 13 (A-Ss 2 12 )e 

H(t)~ 5c 12 s 12 5c\ 2 -9 13 5c 12 s 12 e- iS ^ 

K 9 13 (A-5s 2 12 )e i5 ° p -9 13 5c 12 s 12 e i5 ^ A 

Separating H(t) in independent and dependent parts of 9 i3 yields 

H(t) = H (t) + H', H' = H 1 + H 2 , 



i<5cp\ 



(A10) 



(AH) 



34 



where 



's 2 2 6 + V(t) c 12 s 12 5 



H (t) 



C12S12S 




c 2 12 5 




h(t) 



A 



V 



(A12) 



6 13 (A-s 2 l2 5)e~^\ 

-0 13 c l2 s 12 e- i5cp 5 

\e 13 (A-s 2 12 5)e iS °r -9 13 c 12 s 12 e iS °r5 

'() a 
0b 



A 



a* 6* 



# 9 



V 



137- 



(A13) 
(A14) 



Here the Hamiltonian Hi is of order # 13 , whereas the Hamiltonian if 2 is of order 8f 3 . Note 
that the Hamiltonian H 1 is independent of time t. Furthermore, the time-dependent Hamil- 
tonian Ho(t) is only dependent on the mixing angle 612. 

Inserting Eq. (jAll|) as well as S(t,t ) = So(t,to)Si(t,t ) into Eq. ()A8|) gives 



i Qj^o(Mo)) Si(Mo) + LSb(Mo)j|s , i(*,*o) 
H (t)S (t, t )Si(t, t ) + #iSb(t, t )Si(t, to)- 



(A15) 



Now, assuming that i4rSo(t, to) = H (t)So(t,to) holds implies that we have the equation 
i^Si(t,to) = Hi(t)Si(t,to), where i^i(t) = 5 , _1 (t, t )HiSo(t, to), which can be integrated to 
give the integral equation 

^(Mo) = l-i [ Hi(t')Si(t',t )dt> 

J to 

= i-i / 5 - 1 (^ / ,^o) J fflS , o(t / ,to)5 1 (^ / ,^o)dt , . (A16) 



(A17) 



Thus, from first order perturbation theory we obtain |43j , |4 



S(t, t ) ~ 5 (t, t ) - i5 (t, t ) / S to^Soit', t )dt'. 

'to 



Since we assumed before that i^So(£, to) = Ho(t)So(t,to) holds, we have now to find 
5*o(t, t ). We observe that the 2x2 submatrix in the upper-left corner of H (t) in Eq. (jA12|) . 
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i.e., h(t), is not traceless. Making this submatrix traceless yields 



H (t) = H (t)--txh(t)l 3 

f -Wl 2 -s\ 2 )8+\V{t) 
C12S128 




Cl2Sl2<5 



Hc\ 2 -s\ 2 )5-\V{t) 



2 V-12 









\ 



A-±[V(t)+5]) 



(A18) 



Note that, in general, any term proportional to the identity matrix I3 can be added to or 
subtracted from the Hamiltonian H (t) without affecting the neutrino oscillation probabil- 
ities. In particular, a term such that the 2x2 submatrix h(t) in the upper-left corner of 
Ho(t) becomes traceless [see Eq. (jA18|) ]. Furthermore, note that the new Hamiltonian Ho(t) 
will not be traceless and that the (3, 3)-element of H (t) will, of course, also be changed by 
such a transformation. 

Instead of solving i^So(t, to) = H (t)So(t,to), we have now to solve i^Soft, t ) = 
Hn(t) Sn( t : t n ) + ^t r h(t)Sn(t, to). The solution to this equation, So(t, t ), has the general 



form 



43, 



131 



/ 



So{t, tc 



\ 



V 



(A19) 



a(t,t ) p(t,t ) 
-P*(t,to) a*(t,to) 
f{t,t )J 

where the functions a(t, to) and /3(t,to) describe the two flavor neutrino evolution in the 
(1, 2)-subsector, in which the 2x2 submatrix h of H acts as the Hamiltonian. In the end of 
this appendix, we will derive the analytical expressions for the functions a(t, t ) and f3(t, t ). 
The function f(t,t ) can, however, immediately be determined to be 



/(Mo) 



-ij* A(t'W 



(A20) 



where A(t) = A - \ [V(t) + 5} and 



$(t,t 



A(t')dt' 



to 

A 



to 



1 r A-i[V(t') + 5]K 



V(t')dt'. 



to 



36 



Now, inserting Eqs. (lA"T3j) and (jMgj) into Eq. (lATH) yields 

' a(t,t ) /3(t,t ) -if(t,t )A(t,t ) 

Si I; t ) ~ to) «*(t,t ) -if(t,t )B(t,t ) 

-if{t,t )C(t,t ) -if{t,t )D{t,t ) /(t,t ) 



V 



(A21) 



where 



A(t,t ) = f*(t,t Q ){a[a(t,t )I Q ., t0 (t,t ) + P(t,t Q )I p ., t0 (t,t )] 

+ b[p(t,t )I a)to (t,t ) - a(t,t )If3 )to (t,to)) }, 

B{t,t ) = f*(t,t ){a[a*(t,t Q )Ip*, t0 (t,t ) - P*(t,t )I a *, to (t,t Q )} 

+ b[a*(t,t Q )I a>t0 (t,t ) +/3*(t,t )If),t (t,t )] }, 

C(Mo) = a*r a ^ to (t,t )-b*Il to (t,t ), 

D(t,t Q ) = a*i;* >t0 (t,t ) + b*I* ajt0 (t,t Q ) 



with 



I<p,to{ti ^o) 



<p(t',t )f(t',t )dt', p = a,a*,P,P*. 



(A22) 

(A23) 
(A24) 
(A25) 

(A26) 



to 



Equations (jA22|) and (jA23|) can be further simplified using the following 



£>o(ti,t) — So(ti,to)S (t,to) 
( 



a(t 1; t ) P(t u t ) 
-/3*(h,t ) a* (Mo) 
f{h,t )J 

\ 



a*(t,t ) -P(t,to) 
P*(t,t ) a(t,t ) 
y f*(t,to)J 



a(t u t) P(h,t) 
-P*(t 1} t) a*(t u t) 
\ o f{h,t)J 

Considering Eq. (jA27|) . one immediately finds that 

a(ti,*o)a*(Mo) + P{ti,t )P*{t,t ) = a(h,t), 

-a(ti,t )P(t,t ) + P(ti,t )a(t,t ) = P(h,t), 

-/F(Mo)a*(Mo) + a*(Mo)/3*(Mo) = ~P*{h,t), 

a*(t u t )a(t,t ) + P*(t 1 ,t Q )P(t,t ) = ct{t u t), 

f(t x ,to)r{t,to) = f(h,t). 



(A27) 



(A28) 
(A29) 
(A30) 
(A31) 
(A32) 
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Thus, using Eqs. jggfl - (lA"32j) as well as the identity \ f(t,t )\ 2 = f(t,t )f*(t,t ) = 1, one 
can write Eqs. (jA22j) and (jA23|) as 



A(t,t ) = al a * )t (t,t ) -blp )t (t,t ), 
B(t,t ) = alp*, t (t,t )+bl a>t (t,t ). 



(A33) 
(A34) 



Now, rotating S(t, to) back to the original basis, one finds the evolution operator for neutrinos 
in the flavor basis 

S f (t,t ) = O^ 3 S(t,t )O 23 

1 a C23(3-is 23 fA s 23 (3 -ic2 3 fA 

-C23(3* ~is 2 3fC S 2 2 S 2 3 =(Sf,ab), 

^ S23P* ~ lc 23fC S32 S33 j 



(A35) 



where 



(A36) 
(A37) 
(A38) 
(A39) 

B(t,t ), 



522 = cl 3 a* + s 2 23 f - ls 2 3C 2 3f(B + D), 

523 = -S23C23 O* - /) - i/ (c 2 23 B - s 2 23 D) , 
5*32 = -S23C23 (a* - f) + if (sl 3 B - c\ 3 D) , 
S33 = sj 3 a* + c 2 23 f + is 23 C23/ (B + D) 

with the notation a = a(t,t ), (3 = j3(t,t ), f = f(t,t ), A = A(t,t ), B 
C = C{t,t ), and D = D(t,t ). 

Similarly, replacing the total Hamiltonian for neutrinos (jA4|) with the total Hamilto- 
nian for antineutrinos ()A5j) in the Schrodinger equation (|A3J) . the evolution operator for 
antineutrinos in the flavor basis becomes 

( a C23/3 - is 23 fA -S23P - ic 2 3fA* 

-C23P* ~ IS23JC S 2 2 S 2 3 = (Sf,ab) 

l S23P* — lC23fC S32 S33 



Sf(t, t c 



(A40) 



/ 



where 



5 22 = c 23 a* + s 23 f - is 2 3C 2 3f(B + D), 

523 = -S23C23 (a* - f) - if (c 2 23 B - s 2 23 D) , 
5*32 = -S23C23 (a* - f) + if (sj 3 B - c 2 23 D) , 
£33 = s\ 3 a* + c 2 23 f + 1S23C23J (B + D) 
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(A41) 
(A42) 
(A43) 
(A44) 



with the same type of notation as in the neutrino case. 

We will now derive the general analytical expressions for the functions a(t, to) and /3(t, to). 
In order to perform this derivation, we study the evolution operator in the (1, 2)-subsector, 
which is a separate problem in the rotated basis, and its solution is independent from the 
total three flavor neutrino problem. We assume that the evolution operator in the (1, 2)- 
subsector, S(i t2 )(t, to), satisfies the Schrodinger equation for neutrinos 

^#(1,2) (Mo) = h(t)S (h2) (t,to), (A45) 
where h(t) is the Hamiltonian and it is given by 
s\ 2 5 + V(t) s 12 c 12 5 



hit) 



s 12 c 12 5 c\ 2 5 

-\&-^6+m ^ ) + 1 - [ s + v(t) ] i 2 , 



(A46) 



see Eqs. (IA12)) and (|A18j) . Note that the term proportional to the identity matrix 1 2 in the 
Hamiltonian h(t) does not affect neutrino oscillations, since such a term will only generate 
a phase factor. Thus, we need not consider this term. In addition, note that the same 
term has been subracted from the Hamiltonian H (t) [see Eq. (jA18|) ] for the total three 
flavor neutrino problem. Thus, it also in this case only gives rise to a phase factor in the 
three flavor neutrino evolution operator So(t, to) [see Eq. (IA19J) ]. which does not affect the 
neutrino oscillations. 

The solution to the Schrodinger equation in the (1, 2)-subsector is 

5(i,2) (t, to) = e-< h{t ' )dt ' = e- iH ^ t0 \ (A47) 
where the integrated Hamiltonian, H(t,t ), is given by 

1 / - cos 26 12 5(t - t ) + f! V(t')dt' sin 26 12 5(t - t ) \ 

if (Mo) = - t0 • (A48) 

2 \ sin 26 12 5(t - t ) cos 29 12 S(t - t ) - /* V(?)dtJ 

Since H(t,to) is a 2 x 2 matrix, the solution can be written on the following form 



S(i,2)(Mo) = cos v / -detif(t,t )l 2 



1 



- i— ====== sin y/- det Hit, t )H(t, t ), (A49) 

y/-destH{t,t Q ) 
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where the determinant of H(t,to), det H(t, to), is given by 



det if (Mo 



cos26 12 5{t-t ) - [ V(t')dt 

Jt 



-i 2 



-sin 2 2^ 12 (5 2 (t-to 



+ 



S 2 (t-t ) 2 -2 cos 26 12 5(t- t ) / V(t')dt' 

Jto 

2 ■ 

V{t')dt> ' 



to 



(A50) 



Furthermore, the eigenvalues of H (t, t ) can be found from the characteristic equation 
det (H (t, to) — fil 2 ) = 0, which yields f2 = ±a/— det H(t, to). Note that in vacuum, i.e., 
V(t) = W, it holds that fi 2 |y( t ) =0 = |5 2 (t — t ) 2 = f2 2 ac . Now, if one writes the evolution 
operator S/i )2 )(t, to) as 



£(1,2) (Mo) 



a(t,t ) *o) 
-/3*(Mo) a*(Mo) 



(A51) 



then, using Eq. ()A49|) . one can identify the functions a(t, t ) and j3(t,t Q ). We obtain 

.sinfi 



a(t, t ) = cosf2 + i- 



2Q 



cos20i 2 5(t-t o ) - / V(*) d * 

to 



sin Q 

(3{t,t ) = -i— S m26 12 5{t - t ) 



(A52) 
(A53) 



where again 



n = ±y/- det H(t, t 
, S(t - t 



COS 2^12 — 



1 



V{t')df 



S(t ~ to) Jto 

Similarly, for antineutrinos the functions a(t,t ) and /3(t,t ) become 

.sinO 



sin 2 26 



12- 



(A54) 



a(t, t ) = cosfi + i- 



2Q 



cos28 12 5(t-to) + I V(t')dt' 

to 



— sin Q 

(3(t,t ) = -i_sin2# 12 5(t-t ), 



(A55) 
(A56) 



which we, in principle, obtain by making the replacement V(t) — > —V(t) in the expressions 
for the functions a(t,to) and fi(t,to). Here 



S(t - to) 



cos20 12 + 



1 



S(t-t ) J t0 



V{t')dt' 



+ sin 2 26*12. 



(A57) 
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Note that ft in Eq. (|A54|) and ft in Eq. (|A57|) only differ with respect to the sign in front of 
the integral of the matter potential. Thus, from the expressions for ft and ft we find that 

(A58) 



n 2 = ft 2 + cos26 12 5(t-t ) I V(t')dt'. 

Jt 



Let us now consider some special cases when the relation between ft and ft becomes simpler. 
In the case that 



t — to = 0, one finds ft = ft = 0, which is a trivial and non-interesting case. 

S = 0, we have degenerated neutrino masses mi = m 2 (and negligible solar mass 
squared difference) or extremely high neutrino energy and this leads to ft 2 = ft 2 = 



to 

i§ lr fta v (t') dt ', « = cosft + i§|p j'l V(t')dt' = a*, and (3 = = 0. 



f tn V(t')dt' , which implies that ft = ±ft. Thus, in addition, we have a = cos ft 

inQ ft 
2Q Jto 



cos2# 12 = (e.g., 6 12 = 45°), we have maximal mixing in the (1, 2)-subsector and this 



1 2 



leads to ft 2 = ft 2 = \5 2 {t - t ) 2 + ~ J£ V(t')dt' , which again implies that ft = ±ft. 



a 



In this case, we find a = cos ft - i^fp J t [ V{t')dt', a = cos ft + i^fp j£ V{t')dt' 
and/5 = ^ = -i^5{t-t ). 

It V(t')dt' = 0, one obtains ft 2 = ft 2 = \5 2 (t — t ) 2 , which also implies that ft 
±ft. Furthermore, one has a = a = cos ft + i^fpp cos 29i 2 S(t — t ) and j3 = j3 



• sin Q 

2n 



sm26 12 5(t-t ). 



In addition, if we have close to maximal mixing, i. e., Q\ 2 < 45°, then we can write 812 = f — e, 
where e is a small parameter. Making a series expansion with the parameter e as a small 
expansion parameter, we obtain 

4 



cos 26*12 
ft 



2e - -e 3 + 0{e 5 



± 



ft + ^Sit - t ) 



V{t')dt'e + C(e 2 



to 



(A59) 
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[1] S. R. Coleman and S. L. Glashow, Phys. Lett. B405, 249 (1997), hep-ph/9703240, 
[2] S. R. Coleman and S. L. Glashow, Phys. Rev. D59, 116008 (1999), |hejTph/9812418| 



41 



[3] V. D. Barger, S. Pakvasa, T. J. Weiler, and K. Whisnant, Phys. Rev. Lett. 85, 5055 (2000), 



hep-ph/0005197 



[4] 
[5] 

[6] 

[7] 



[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17; 

[18 
[19 
[20 

[21 

[22 
[23 

[24 

[25 



H. Murayama and T. Yanagida, Phys. Lett. B520, 263 (2001), hep-ph/0010178. 

M. C. Bahuls, G. Barenboim, and J. Bernabeu, Phys. Lett. B513, 391 (2001), 

|hep-ph/0102"l84l 

G. Barenboim, L. Borissov, J. Lykken, and A. Y. Smirnov, J. High Energy Phys. 10, 001 
(2002), |hep"-ph/0108199| 

G. Barenboim, L. Borissov, and J. Lykken, Phys. Lett. B534, 106 (2002), piep^ph/ 0201080, 
G. Barenboim, J. F. Beacom, L. Borissov, and B. Kayser, Phys. Lett. B537, 227 (2002), 
|hep-ph/0203"26l 



G. Barenboim and J. Lykken, Phys. Lett. B554, 73 (2003), |heiTph/0210411| 



G. Barenboim, L. Borissov, and J. Lykken, hep-ph/0212116 
S. Pakvasa, |hep-ph/0110175l 



Z.-z. Xing, J. Phys. G28, B7 (2002), hep-ph/0112120 



S. Skadhauge, Nucl. Phys. B639, 281 (2002), [hepjjh/0112189 

S. M. Bilenky et al, Phys. Rev. D65, 073024 (2002), jhep-ph/0112226 

T. Ohlsson, hep-ph/0209150, 



J. N. Bahcall, V. Barger, and D. Marfatia, Phys. Lett. B534, 120 (2002), hep-ph/0201211 
A. Strumia, Phys. Lett. B539, 91 (2002), [hep-ph/0201134[ 



I. Mocioiu and M. Pospelov, Phys. Lett. B534, 114 (2002), |hef>ph7 0202160 . 



A. De Gouvea, Phys. Rev. D66, 076005 (2002), hep-ph/0204077 



LSND Collaboration, C. Athanassopoulos et al, Phys. Rev. Lett. 77, 3082 (1996), 
|nucl-ex/ 9605003 

LSND Collaboration, C. Athanassopoulos et al, Phys. Rev. Lett. 81, 1774 (1998), 
|nucl-ex/ 9709006 

LSND Collaboration, A. Aguilar et al, Phys. Rev. D64, 112007 (2001), pIep^ex/0104049 



Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998), 
|hep-ex/ 9807003. 

Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 82, 2644 (1999), 
|hep-ex/9812014l 

Super-Kamiokande Collaboration, S. Fukuda et al, Phys. Rev. Lett. 85, 3999 (2000), 



42 



[26 
[27 



[29 
[30 

[31 

[32 

[33 

[34 

[35 

[36 

[3r 



|hep-ex/ 0009001 . 

M. Shiozawa, talk given at the XXth International Conference on Neutrino Physics & As- 
trophysics (Neutrino 2002), Munich, Germany, 2002. 

Super-Kamiokande Collaboration, S. Fukuda et at, Phys. Rev. Lett. 86, 5651 (2001), 
|hep-ex/0103032" 



Super-Kamiokande Collaboration, S. Fukuda et al., Phys. Rev. Lett. 86, 5656 (2001), 
|hep-ex/0103"033l 

Super-Kamiokande Collaboration, M. B. Smy, hep-ex/0106064 

Super-Kamiokande Collaboration, S. Fukuda et al, Phys. Lett. B539, 179 (2002), 
hep-ex/0205075 



M. Smy, talk given at the XXth International Conference on Neutrino Physics & Astrophysics 
(Neutrino 2002), Munich, Germany, 2002. 

SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 87, 071301 (2001), 
|nucl-ex/0106015 



SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002), 
|nucl-ex/ 0204008 

SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 89, 011302 (2002), 
|nucl-ex/020 4009 

A. Hallin, talk given at the XXth International Conference on Neutrino Physics & Astro- 
physics (Neutrino 2002), Munich, Germany, 2002. 
MiniBooNE Collaboration, http://www-boone.fnal.gov/ 



J. N. Bahcall, M. C. Gonzalez-Garcia, and C. Peha-Garay, J. High Energy Phys. 07, 054 



(2002), hep-ph/0204314 



[38] G. Liiders, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 28, 5 (1954). 

[39] W. Pauli, in Niels Bohr and the Development of Physics, edited by W. Pauli, L. Rosenfeld, 

and V. Weisskopf (Pergamon, London, 1955). 
[40] J. S. Bell, Proc. R. Soc. London A231, 479 (1955). 

[41] KamLAND Collaboration, K. Eguchi et al., Phys. Rev. Lett. 90, 021802 (2003), 
|hep-ex/0212021| 

[42] J. N. Bahcall, M. C. Gonzalez-Garcia, and C. Peha-Garay, J. High Energy Phys. 02, 009 
(2003), |hep-ph/0212147 



43 



[43] E. K. Akhmedov, P. Huber, M. Lindner, and T. Ohlsson, Nucl. Phys. B608, 394 (2001), 
|hep-ph/0105029[ 



[44; 
[45 

[46 
[47; 



[49 

[50; 

[51 
[52 

[53 

[54 
[55 
[56 

[57; 
[58 



[60 
[61 

[62 

[63 
[64 



J. Bernabeu, hep-ph/9904474 

J. Bernabeu, S. Palomares-Ruiz, A. Perez, and S. T. Petcov, Phys. Lett. B531, 90 (2002), 
|hep-ph/0110"07T| 

J. Bernabeu and S. Palomares-Ruiz, hep-ph/0112002 



J. Bernabeu and S. Palomares-Ruiz, Nucl. Phys. Proc. Suppl. 110, 339 (2002), 
|hep-ph/020i090l 



J. Arafune, M. Koike, and J. Sato, Phys. Rev. D56, 3093 (1997), hep-ph/9703351 



H. Minakata and H. Nunokawa, Phys. Rev. D57, 4403 (1998), piep^ph /9705208 
H. Minakata and H. Nunokawa, Phys. Lett. B413, 369 (1997), piep^ph /9706281 . 
K. R. Schubert, |hep"-ph/9902215 



K. Dick, M. Freund, M. Lindner, and A. Romanino, Nucl. Phys. B562, 29 (1999), 
|hep-ph/9903"308 



A. Donini, M. B. Gavela, P. Hernandez, and S. Rigolin, Nucl. Phys. B574, 23 (2000), 
|hep-ph/ 9909254 



A. Cervera et al, Nucl. Phys. B579, 17 (2000), hep-ph/0002108 ; B593, 731(E) (2001). 

H. Minakata and H. Nunokawa, Phys. Lett. B495, 369 (2000), |h"ei>ph/0004114l 

V. D. Barger, S. Geer, R. Raja, and K. Whisnant, Phys. Rev. D63, 033002 (2001), 

|hep-ph/000718T 



P. M. Fishbane and P. Kaus, Phys. Lett. B506, 275 (2001), hep-ph/0012088 



M. C. Gonzalez-Garcia, Y. Grossman, A. Gusso, and Y. Nir, Phys. Rev. D64, 096006 (2001), 
|hep-ph/0105"T59"l 

H. Minakata and H. Nunokawa, J. High Energy Phys. 10, 001 (2001), |hei>ph/0108085 



H. Minakata, H. Nunokawa, and S. Parke, Phys. Lett. B537, 249 (2002), hep-ph/020417T| 
K. Kimura, A. Takamura, and H. Yokomakura, Phys. Rev. D66, 073005 (2002), 
|hep-ph/0205295l 

H. Yokomakura, K. Kimura, and A. Takamura, Phys. Lett. B544, 286 (2002), 
|hep-ph/0207174l 

T. Ota and J. Sato, Phys. Rev. D67, 053003 (2003), hep-ph/0211095l 
T. K. Kuo and J. Pantaleone, Phys. Lett. B198, 406 (1987). 



44 



[65] P. I. Krastev and S. T. Petcov, Phys. Lett. B205, 84 (1988). 

[66] S. Toshev, Phys. Lett. B226, 335 (1989). 

[67] S. Toshev, Mod. Phys. Lett. A6, 455 (1991). 

[68] J. Arafune and J. Sato, Phys. Rev. D55, 1653 (1997), |heiTph/9607437| 

[69] M. Koike and J. Sato, |hep-ph/9707"203l 

[70] S. M. Bilenky, C. Giunti, and W. Grimus, Phys. Rev. D58, 033001 (1998), pIep^ph/9712537 



[71] V. D. Barger, Y.-B. Dai, K. Whisnant, and B.-L. Young, Phys. Rev. D59, 113010 (1999), 
|hep-ph/9901388[ 



[72] M. Koike and J. Sato, Phys. Rev. D61, 073012 (2000), hep-ph/9909469. 

[73] O. Yasuda, Acta Phys. Polon. B30, 3089 (1999), ||hep-ph/9910428[ 

[74] J. Sato, Nucl. Instrum. Meth. A451, 36 (2000), hep-ph/ 99 104421 

[75] M. Koike and J. Sato, Phys. Rev. D62, 073006 (2000), plep^"ph/9911258| 

[76] P. F. Harrison and W. G. Scott, Phys. Lett. B476, 349 (2000), |hejTph/9912435[ 

[77] V. A. Naumov, Sov. Phys. JETP 74, 1 (1992), zh. Eksp. Teor. Fiz. 101 (1992) 3. 

[78] V. A. Naumov, Int. J. Mod. Phys. Dl, 379 (1992). 

[79] H. Yokomakura, K. Kimura, and A. Takamura, Phys. Lett. B496, 175 (2000), 
hep-ph/0009141 



[80] S. J. Parke and T. J. Weiler, Phys. Lett. B501, 106 (2001), hep-ph/0011247. 



[81] M. Koike, T. Ota, and J. Sato, Phys. Rev. D65, 053015 (2002), hep-ph/0011387 



[82] T. Miura, E. Takasugi, Y. Kuno, and M. Yoshimura, Phys. Rev. D64, 013002 (2001), 
|hep-ph/0102111 



T. Miura, T. Shindou, E. Takasugi, and M. Yoshimura, Phys. Rev. D64, 073017 (2001), 
|hep-ph/0106086| 
[84] A. Rubbia,|hep^ph70106088 



[85] Z.-z. Xing, Phys. Rev. D64, 093013 (2001), [hep-ph/0107005 



[86] T. Ota, J. Sato, and Y. Kuno, Phys. Lett. B520, 289 (2001), hep-ph/0107007 
[87] T. Ohlsson, hep-ph/0108048. 



[88] A. Bueno, M. Campanelli, S. Navas-Concha, and A. Rubbia, Nucl. Phys. B631, 239 (2002), 
|hep-ph/0112297 



J. Sato, eConf C010630, E105 (2001). 
[90] E. K. Akhmedov, |hep-ph/ 0207342} 



45 



[91] H. Minakata, H. Nunokawa, and S. Parke, Phys. Rev. D66, 093012 (2002), hep-ph/0208163. 
[92] P. Ruber, J. Phys. G29, 1853 (2003), plep^"ph/0210140| 

[93] C. N. Leung and Y. Y. Y. Wong, Phys. Rev. D67, 056005 (2003), |hejrph/0301211[ 
[94] S. M. Bilenky and B. Pontecorvo, Phys. Rept. 41, 225 (1978). 

[95] S. M. Bilenky and S. T. Petcov, Rev. Mod. Phys. 59, 671 (1987), 61 (1989) 169(E). 
[96] M. Lindner, hep-ph/0209083 and hep-ph/0210377, and references therein. 
[97] H. Minakata, private communication. 
[98] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985). 
[99] C. Jarlskog, Z. Phys. C29, 491 (1985). 
[100] J. N. Bahcall, M. C. Gonzalez-Garcia, and C. Peha-Garay, J. High Energy Phys. 08, 014 

(2001) , |hep-ph/010"6258l 

[101] J. N. Bahcall, M. C. Gonzalez-Garcia, and C. Peha-Garay, J. High Energy Phys. 04, 007 

(2002) ,|hep-ph/01lTT50l 



[102] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B420, 397 (1998), |hep^ex7971 1002 



[103] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B466, 415 (1999), hep-ex/ 9907037, 
[104] CHOOZ Collaboration, M. Apollonio et al, Eur. Phys. J. C27, 331 (2003), hep-ex/0301017l 
[105] A. M. Dziewonski and D. L. Anderson, Phys. Earth Planet. Interiors 25, 297 (1981). 



[106] M. Freund and T. Ohlsson, Mod. Phys. Lett. A15, 867 (2000), hep-ph/9909501 , and refer- 
ences therein. 

[107] BNL Neutrino Working Group Collaboration, D. Beavis et al, hep-ex/0205040. 



[108] W. J. Marciano, hep-ph/0108181 



[109] M. V. Diwan et al, Phys. Rev. D68, 012002 (2003), |hep^ph /0303081 

[110] BooNE Collaboration, E. D. Zimmerman, eConf C0209101, TH05 (2002), hep-ex/0211039l 

[111] BooNE Collaboration, A. O. Bazarko, |hep-ex/0210 020 

[112] ICARUS Collaboration, A. Rubbia et al, cERN-SPSLC-96-58. 

[113] ICARUS Collaboration, F. Arneodo et al, hep-ex/ 0103008 

[114] OPERA and ICARUS Collaboration, D. Duchesneau, eConf C0209101, TH09 (2002), 
|hep-ex/0209082[ 

[115] JHF-Kamioka Collaboration, Y. Itow et al, |hep-ex/ 0106019 , 



[116] K2K Collaboration, S. H. Ahn et al, Phys. Lett. B511, 178 (2001), [hej>ex/0103001 . 



[117] K2K Collaboration, M. H. Ahn et al, Phys. Rev. Lett. 90, 041801 (2003), |hep^ex70 212007 

46 



[118] MINOS Collaboration, S. G. Wojcicki, Nucl. Phys. Proc. Suppl. 91, 216 (2001). 

[119] MINOS Collaboration, V. Paolone, Nucl. Phys. Proc. Suppl. 100, 197 (2001). 

[120] MINOS Collaboration, M. V. Diwan, eConf C0209101, TH08 (2002), |hep^ex/0211026| 



[121] NuMI Collaboration, D. Ayres et al, hep-ex/0210005 



[122] NuTeV Collaboration, S. Avvakumov et al, Phys. Rev. Lett. 89, 011804 (2002), 
|hep-ex/0203018[ 

[123] OPERA Collaboration, H. Pessard, Phys. Scripta T93, 59 (2001). 

[124] Palo Verde Collaboration, F. Boehm et al, Prog. Part. Nucl. Phys. 40, 253 (1998). 



[125] Palo Verde Collaboration, F. Boehm et al, Phys. Rev. Lett. 84, 3764 (2000), hep-ex/9912050 
[126] Palo Verde Collaboration, F. Boehm et al, Phys. Rev. D62, 072002 (2000), lhep-ex/ 0003022 
[127] Palo Verde Collaboration, F. Boehm et al, Phys. Rev. D64, 112001 (2001), hep-ex/0107009| 
[128] T. Ohlsson and H. Snellman, J. Math. Phys. 41, 2768 (2000), |hep^ph/9910546t D42, 2345(E) 
(2001). 

[129] T. Ohlsson and H. Snellman, Phys. Lett. B474, 153 (2000), |heiTph/9912295| B480, 419(E) 
(2000). 

[130] T. Ohlsson, Phys. Scripta T93, 18 (2001). 

[131] E. K. Akhmedov, Nucl. Phys. B538, 25 (1999), |hep^ph/9805272 
[132] E. K. Akhmedov, |hep-ph/0001264l 

[133] Particle Data Group, K. Hagiwara et al., Phys. Rev. D66, 010001 (2002), 
|http://pdg. lbl.gov/ 



[134] For three neutrino flavors, we use the stan dard parameterization of the unitary leptonic 



mixing matrix U = £/(6>i2, 6>i 3 , # 2 3> #Cp) 



1331 ] , where Sj,- = sin 9™ and = cos 9™ . Here 9\2 , 



#13, and #23 are the leptonic mixing angles and 5qp is the leptonic CP violation phase. 
[135] The mantle-core-mantle step-function approximation of the Earth matter density profile 
consists of three constant matter density layers (mantle, core, and mantle) with p man tie = 
4.5g/cm 3 and p core = 11.5g/cm 3 . 



[136] We will use the standard parameterization of the leptonic mixing matrix 



133] 



47 



